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Abstract 


Within the framework of the theory of the column and row determi¬ 
nants, we obtain new determinantal representations of the W-weighted 
Drazin inverse over the quaternion skew field. We give determinantal 
representations of the W-weighted Drazin inverse by using previously in¬ 
troduced determinantal representations of the Drazin inverse, the Moore- 
Penrose inverse, and the limit representations of the W-weighted Drazin 
inverse in some special case. 

1 Introduction 

Let R and C be the real and complex number fields, respectively. Throughout 
the paper, we denote the set of all m x n matrices over the quaternion algebra 

H = {oo + aii + a 2 j + ask \ f = —1, oq, oi, 02,03 S R} 

by and by the set of all m x n matrices over HU with a rank r. 

Let M (n, H) be the ring of n x n quaternion matrices and I be the identity 

matrix with the appropriate size. For A G we denote by A*, rank A the 

conjugate transpose (Hermitian adjoint) matrix and the rank of A. The matrix 
A = (a^) G is Hermitian if A* = A. 

The definitions of the generalized inverse matrices may be extended to 
quaternion matrices. 

The Moore-Penrose inverse of A G denoted by A^, is the unique 

matrix X G satisfying the following equations. 


AXA = A; 
XAX = X; 
(AX)* = AX; 
(XA)* = XA. 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


For A G with k = Ind A the smallest positive number such that rank A^+^ 

rankA^ the Drazin inverse of A, denoted by , is defined to be the unique 
matrix X that satisfying (1.2) and the following equations. 


AX = AX; 

A'=+iX = A^ 


( 5 ) 

( 6 ) 
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In particular, when IndA. = 1, then the matrix X is called the group inverse 
and is denoted by X = A®. 

If IndA = 0, then A is nonsingular, and A^ = A^ = A~^. 

Cline and Greville [T] extended the Drazin inverse of square matrix to rect¬ 
angular matrix, which can be generalized to the quaternion algebra as follows. 
For A € and W S the W-weighted Drazin inverse of A with 

respect to W, denoted by A^ vVi is the unique solution to equations. 


(AW)'=+iXW = (AW)'=; (7) 

XWAWX = X; (8) 

AWX = XWA, (9) 

where k = max{/n(i(AW),/nd(WA)}. It is denoted by X = Ad^w- The 
properties of the complex W-weighted Drazin inverse can be found in [IHB]. 
These properties can be generalized to H. If A G with respect to W G 

Ijnxm ^ niax{/nd(AW),/nd(WA)}, then 

Ad,w = A((WA)^))' = ((AW)^))'A, (10) 

Ad,wW = (WA)^), WAd.w = (AW)^. (11) 


The problem of determinantal representation of generalized inverse matrices 
only recently begun to be decided through the theory of the column-row de¬ 
terminants introduced in [HIS]. The theory of row and column determinants 
develops the classical approach to a definition of a determinant, as alternating 
sum of products of the entries of matrix but with a predetermined order of 
factors in each terms of the determinant. A determinant of a quadratic matrix 
with noncommutative elements is often called the noncommutative determi¬ 
nant. Unlike other known noncommutative determinants such as determinants 
of Dieudonne Study [To], Moore mm. Chen [13], quasideterminants of 
Gelfand-Retakh [14] , the double determinant built on the theory of the column- 
row determinants has properties similar to a usual determinant, in particular it 
can be expand along arbitrary rows and columns. This property is necessary 
for determinantal representations of an inverse and generalized inverse matri¬ 
ces. Determinantal representations of the Moore-Penrose inverse and the Drazin 

inverse over the quaternion skew-field have been obtained in [151116] and m. 

( 2 ) 

respectively. Determinantal representations of an outer inverse A^g is intro¬ 
duced in [121 [IS] using the column-row determinants as well. Recall that an 
outer inverse of a matrix A over complex field with prescribed range space T 
and null space S' is a solution of (1.2) with restrictions, 

7^(X) = T, A/'(X) = S. 

Within the framework of the theory of the column-row determinants Song [20] 
also gave a determinantal representation W-weighted Drazin inverse over the 
quaternion skew-field using a characterization of the W-weighted Drazin inverse 
by an outer inverse A^ g. But in obtaining of this determinantal representation 
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is used auxiliary matrices which different from A or its powers. In this paper 
we obtain determinantal representations of the W-weighted Drazin inverse of 
A € with respect to W € by using only their entries. 

The paper is organized as follows. We start with some basic concepts and 
results from the theory of the row and column determinants and give the de¬ 
terminantal representations of the inverse, the Moore-Penrose inverse, and the 
Drazin inverse over the quaternion skew field in Section 2. In Section 3, we 
obtain determinantal representations of the W-weighted Drazin inverse by us¬ 
ing introduced above determinantal representations of the Drazin inverse, the 
Moore-Penrose inverse, and the limit representations of the W-weighted Drazin 
inverse in some special case. In Section 4, we show a numerical example to 
illustrate the main result. 


2 Elements of the theory of the column and row 
determinants 

For a quadratic matrix A = (aij) € M (n, H) we define n row determinants and 
n column determinants as follows. Suppose is the symmetric group on the 
set In = { 1 , • ■ • ,n}. 

Definition 2.1 The ith row determinant of A = (atj) € M(n,]HI) is defined for 
all i = l,n by putting 

rdetiA = (~ 1 ) O,***,! + l + > 

(T^Sn 

^ — {i ik-^ik-i + l • ■ - ^fei-t-Zi) (^fc2^fc2-t-l ■ ■ •^fc2-t-/2) ■ ■ * i.'l’kr-'ikr + l ■ • • ’Ik^ + lr) ? 

with conditions ik^ < iks, < ■ ■ ■ < ikr o,nd ikt < ikt+s for t = 2,r and s = 1 , It- 

Definition 2.2 The jth column determinant of A = (atj) G M(n,]HI) is defined 
for all j = 1, n by putting 

cdetj A = (~1) + jki+li ' ' ' ^jki+llki ^jkiT 

res„ 

~ {jkr+lr ■ ■ ■ jkr + ljkr) ■ ■ ■ (.jk-i+l^ ’ ' ■ jk^ + ljk^) OfelS-il • • • Jfcl-|-1 Jfcl J) ) 

with conditions, jk^ < jk^ < ■ ■ ■ < jk,. and jkt < jkt+s for t = 2,r and s = l,lt- 

Suppose A*^ denotes the submatrix of A obtained by deleting both the tth 
row and the jth column. Let a. j be the jth column and be the ith row of 
A. Suppose A j (b) denotes the matrix obtained from A by replacing its jth 
column with the column b, and A^. (b) denotes the matrix obtained from A by 
replacing its ith row with the row b. 

The following theorem has a key value in the theory of the column and row 
determinants. 
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Theorem 2.1 // A = (oy) £ M (n, H) is Hermitian, then rdetiA = • • • = 

rdet„A = cdetiA = • • • = cdet„A e K. 

Since all column and row determinants of a Hermitian matrix over IHI are equal, 
we can define the determinant of a Hermitian matrix A S M (n, H). By defini¬ 
tion, we put det A := rdeti A = cdeti A, for alH = 1, n. 

The determinant of a Hermitian matrix has properties similar to a usual 
determinant. They are completely explored in by its row and column 

determinants. They can be summarized by the following theorems. 

Theorem 2.2 If the ith row of a Hermitian matrix A € M (n, H) is replaced 
with a left linear combination of its other rows, i.e. a^, = -|- ... -I- Cfca^j.., 

where c; G H for all I = l,k and {i, ii} C In, then 

rdeti Ai. (cia^j. -b ... -b Cfca^j^.) = cdet^ A*, (ciai^. -b ... -b Cfca^j^.) = 0. 

Theorem 2.3 If the jth column of a Hermitian matrix A G M (n, H) is replaced 
with a right linear combination of its other columns, i. e. a,j = aj^ ci -b ... -b 
a.jfcCfe, where c; S H for all I = l,k and {j,ji} C J„, then 

cdetj A j (a -b ... -b a j^Ck) — rdet^ A ^ (a -b ... -b a j,^C}f) — 0. 


The determinant of a Hermitian matrix also has a property of expansion along 
arbitrary rows and columns using row and column determinants of submatri¬ 
ces. So, we were able to get determinantal representations of an inverse and 
generalized inverse matrices as follows. 

Theorem 2.4 ^ If for a Hermitian matrix A G M (n,]HI), 

det A ^ 0, 


then there exist a unique right inverse matrix (RA.) ^ and a unique left inverse 
matrix (LA)~^ of a nonsingular A, where (i?A)~^ = (LA)~^ =: A~^, and the 
right and left inverse matrices possess the following determinantal representa¬ 
tions 


( 12 ) 


Rij = 


j=i i=i 

-rdetj A** (a, i), i^ j, 
rdetfeA”, i=j, 




/ Rii 

i?21 

^ral\ 

(i?A)-' = 

1 

Ri2 

i?22 

Rn2 

det A 

\Rln 

R2n 

^nn J 




/ Lii 

L 21 

Lnl \ 

(LA)-i = 

1 

Li2 

L 22 

Ln2 

det A 

\Lln 

L2n 

Lnn j 


n 

^ij ' ^ij ~ 

n 




Lij — 


( 13 ) 


-cdeti A((' (aj,). 
cdetfc A^ ^ , 


i=j, 
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and A®* (a, i) is obtained from A by both replacing the jth column with the 
ith column and deleting the ith row and column, Aff {blj ,) is obtained by both 
replacing the ith row with the jth row and deleting the jth row and column, 
respectively. In = {1, ■ ■ ■ ,n}, k = min{/„ \ {i}} for all i,j = l,n. 

We shall use the following notations. Let a := {ai,..., at} C { 1 ,..., m} and 
/3 := {/3i,... ,13k} Q {3, ■ ■ ■ ,n} be subsets of the order 1 < fc < min{m,n}. By 
A^ denote the submatrix of A determined by the rows indexed by a and the 
columns indexed by f3. Then A“ denotes the principal submatrix determined by 
the rows and columns indexed by a. If A S M (n,H) is Hermitian, then by |A“| 
denote the corresponding principal minor of det A. For 1 < k < n, the collection 
of strictly increasing sequences of k integers chosen from { 1 ,... ,n} is denoted 
by := {a : a = (ai,..., ak), 3 < ai < ... < ak < n}. For fixed i G a and 
j G fj, let It, ■— } ey I O. G Lr,m^ i G Cy}, Jr,n }j } 13 '• 13 G Lr,n^ j G f3}• 

Denote by a*- and a* the jth column and the ith row of A* and by a*-™^ 

and a^”^^ the jth column and the ith row of A™, respectively. 

The following theorem give determinantal representations of the Moore- 
Penrose inverse over the quaternion skew field H. 

Theorem 2.5 U5f If A G then the Moore-Penrose inverse A+ = 

(aj) G possess the following determinantal representations: 


<■ = 


E cdeD ((A*A) . (a*-)) ^ 

Jr, n.{^} 


E 

3ej, 


(A* A) 


(14) 


or 

Y} rdetj ((AA*)j, (a* )) “ 

E l(AA*)g| • ^ ^ 

CX-^lr, m 

for all i = l,n, j = 1 , m. 

Proposition 2.1 \21} If Ind{A) = k, then A^ = A^(A^^+^)+A^. 

Using the determinantal representations of the Moore-Penrose inverse o and 
di), and Proposition 12.11 we have obtained the following determinantal repre¬ 
sentations of the Drazin inverse for an arbitrary square matrix over H. Denote 
by a.s and a*, the sth column of (A^^+^)*A^ =: A = (a^j G and the ith 

row of A^(A^^+^)* =: A = {hij) G respectively, for all s,t= l,n. 

Theorem 2.6 fl7| / //A G M(n,H) with IndA = k and rank A^+^ = rankA^ = 
r, then the Drazin inverse A^ possess the determinantal representations 


n 


E E cdett 

t=l /3GJr,n{t} 


((Aa+.)-(A»«),(a,))f 


E 

/3eJr, n 


(A2fc+i)* (A2fe+i) ^ 


(16) 
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and 


sf E rdets 

s=i \aG/r,„{4 / .s / y 

E |(A2fe+i(A2fe+i)*) g| 

O'^Ir, n 


(17) 


In the special case, when A S M (n, H) is Hermitian, we can obtain simpler 
determinantal representations of the Drazin inverse. 


Theorem 2.7 [T^IfA £ M (n, H) is Hermitian with IndA = k and rank A^+^ 
rank A^ = r, then the Drazin inverse A^ = (a^) £ possess the following 

determinantal representations: 


„D _ 


E cdet, ((A'^+i) I 

n {^} 


E 

^^Jr, 1 


(Afe+i) 


i\P 


(18) 


or 




E rdetj ((A'=+i)j . (afM “ 

eir,n{j} ^ __ 

E l(A'=+i)g| 

Ct^Ir, n 


(19) 


Note that the determinantal rank of A G (pynxn obtained as the 

largest order of a non-zero principle minor in the Hermitian matrices A* A or 

AA*. 

We shall also need the following facts about the eigenvalues of quaternion 
matrices. Due to the noncommutativity of quaternions, there are two types of 
eigenvalues. A quaternion A is said to be a right eigenvalue of A G M (n, H) 
if A ■ X = X • A for some nonzero quaternion column-vector x with quaternion 
components. Similarly A is a left eigenvalue if A ■ x = A • x for some nonzero 
quaternion column-vector x with quaternion components. The theory on the left 
eigenvalues of quaternion matrices has been investigated in particular in [221124j . 
The theory on the right eigenvalues of quaternion matrices is more developed. 
In particular we note [25H30]. 

Proposition 2.2 \29f Let A G M (u, H) is Hermitian. Then A has exactly n 
real right eigenvalues. 

Right and left eigenvalues are in general unrelated [3T], but it is not for Her¬ 
mitian matrices. Suppose A G M (n, H) is Hermitian and A G M is its right 
eigenvalue, then A-x = x- A = A- x. This means that all right eigenvalues 
of a Hermitian matrix are its left eigenvalues as well. For real left eigenvalues, 
A G K, the matrix AI — A is Hermitian. 

Definition 2.3 If t £ M., then for a Hermitian matrix A £ M(n,]HI) the poly¬ 
nomial PA {t) = det (tl — A) is said to he the characteristic polynomial of A. 
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The roots of the characteristic polynomial of a Hermitian matrix are its real left 
eigenvalues, which are its right eigenvalues as well. We can prove the following 
theorem by analogy to the commutative case (see, e.g. [5^1. 

Theorem 2.8 If A G M(n,]HI) is Hermitian, then p\ (t) = t" — dit^~^ + 
d 2 t^~^ — ... + (—1)"' dn, where dk is the sum of principle minors of A of order 
rk, 1 < k < n, and = det A. 


3 Determinantal representations of the W-weigh- 
ted Drazin inverse for an arbitrary matrix 


Determinantal representations W-weighted Drazin inverse of complex matrices 
have been received by full-rank factorization in [33) and by a limit representation 
in [34] . For an arbitrary matrix over the field of complex numbers, A G 
we denote by 7^(A) the range of A and by A/’(A) the null space of A. For an 
arbitrary matrix over the quaternion skew field, A G we denote by 

TZr{A) = {y G : y = Ax, x G H”}, the column right space of A, 

A/’r(A) = {y £ H" : Ax = 0}, the right null space of A, 

IZi{A) = {y £ H" : y = xA, x £ H™}, the column left space of A, 

Afr{A) = {y £ H™ : xA = 0}, the left null space of A. 

Through the theory of the column-row determinants, a determinantal repre¬ 
sentation W-weighted Drazin inverse over the quaternion skew-field for the first 
time has been obtained in [5D] by the following theorem. 

Theorem 3.1 Let A G W £ withk = ma,x{Ind{AW), IndfWA)} 

and and rank(AW)^ = s. Suppose that B £ and C* £ 

of full ranks and 

TZr{B) = Nr ((WA)'=) , Ar,(C) = Ur {{AWf) , 

7^^(C) = N ((AW)'=), A/((B) = -Ri ((WA)'=) . 

Denote 


M = 


WAW 

C 


B 

0 


Then the W-weighted Drazin inverse A^.w = {o)ij G has the following 

determinantal representations: 


E m+n—s r 

fe=i Lkim^ 


Qij — 


kj 


det M*M 


,i = j = 


( 20 ) 


or 


E m-\-n—s * 7 -, 

k=l ^^kR3k 

det MM* 


= 1,TO, j = l,n, 


( 21 ) 


where Lij are the left (ij)-th cofactor of M*M and Rij are the right (ij)-th 
cofactor of MM*, respectively, for all i,j = l,m n — s. 
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As can be seen, the auxiliary matrices B and C have been used in the deter- 
minantal representations (Uni) and In this paper we escape it. Below we 

give determinantal representations of the W-weighted Drazin inverse of an ar¬ 
bitrary matrix A G with with respect to the matrix W G by using 

the determinantal representations of the Drazin inverse, of the Moore-Penrose 
inverse, and the limit representation of the W-weighted Drazin inverse in some 
particular case. 


3.1 Determinantal representations of the W-weighted 

Drazin inverse by using determinantal representations 
of the Drazin inverse 

Denote WA U = {uij) G and AW =: V = G Due to 

Theorem 12.61 we denote an entry of the Drazin inverse by 


E E cdeD ((U2^+i)* ^ (u.,)) ^ 


D,1 deJr.nft} 


% = 


E (U2fc+l)*(U2fc+l) ^ 

3eJr,n 


( 22 ) 


or 


D,2 

= 


t( E rdet.((u2'=+i(U2fc+i)*) (u,,))“) 

S = 1 \aGlr,r.{s} / .S J J 


,(fe) 

‘'S3 


X; |(U2fc+l (U2'=+l)*) 

aS/r n 


(23) 


where u.g and Ut. the sth column of (u2fc+i)*U^ =: U = (wy) G and the 

rth row of U^(U^*“''^)* =: U = {uij) G H"^", respectively for all s,t = l,n, r = 
rankU^+^ = rankU^. Then by (flUl) we can obtain the following determinantal 
representations of A^.w = 

n 

(24) 

9=1 


where (u^)^^^ = E '^qp^p 3 ^ ^r all l,f= 1,2, and ^nd 

p=i 

from (1251) . Similarly using V = (vy) G we have the following determi¬ 

nantal representations of A^.w, 

m 

(25) 

9=1 
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The first factor is one of the four possible equations 

_ p=i 

alH, / = 1, 2, and an entry of the Drazin inverse is denoting by 


E 


„(fc) 


D,1 /3eJr, m{t} 

= - 


E cdet*((V2^+i)*(V2'=+i) ^(v,)) I 


(v 2 '=+i) (V 2 fc+i) g 

PeJr-.m 


( 26 ) 


or 


D,2 

= 


e( E rdet.((v2'=+i(V^'=+^)*) (v..))“^ 

S=1 \aG/r,m{4 y .s y J 


,(fc) 

Sj 


I (^V2fc+1 (V2*:+l)*) 

OL^lr m 


(27) 


where v.s and v*. the sth column of (V^^+^)*V^ =: V = {vij) G and the 

tth row of V^(V^*^+^)* =: V = {vij) G respectively for all s^t = l,m, 

r = rankV'=+i = rankV'=. 


3.2 Determinantal representations of the W-weighted 

Drazin inverse by using determinantal representations 
of the Moore-Penrose inverse 


Consider the general algebraic structures (GAS) of the matrices A G 
w g A+ e W+ e and Ad,w G with k = 

max{/nd(AW),/n(i(WA)} (e.g., [MS]). Let exist L G and Q € 

such that 


A = L 


All 0 

0 A 22 


Q \ W = Q 


Wii 0 
0 W 22 


Lb 


Then 


O’ 

II 

+ 

■ AE 0 ■ 


W+ = L 

■ WE 

0 ■ 

0 0 



0 

0 


^d,W 


(Wi 


'-11' 

0 


Q 


where Lj Q; All , Wii are non-singular matrices, and A 22 ^^A 22 ,'VV 22 A 22 are 
nilpotent matrices. The follow theorem due to [5] can be expanded to H. 


Theorem 3.2 Let A G and'W G such t/iat A 22 W 22 and'W 22 A 22 

are nilpotent matrices of index k in GAS form. Then the weighted Drazin inverse 
of A with respect to W can he written as matrix expression involving the Moore- 
Penrose inverse, 


Ad,w = {(AW)'= [(AW)^'=+i]^ (AW)'=| W+, (28) 


where k = max{/nd(AW),/nd(WA)}. 
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Similarly can be obtained the following theorem. 

Theorem 3.3 Let A G andW G such t/iat A 22 W 22 and "W 22-1^22 

are nilpotent matrices of index k in GAS form. Then the W-weighted Drazin 
inverse of A with respect to'W can be written as the following matrix expression, 

Ad,w = W+ I (WA)'= [(WA)2 '=+i] + (WA)'= I, (29) 

where k = max{/nc?(AW),/n(i(WA)}. 

Proof. Since W 22 A 22 is a nilpotent matrix of index k, then due to the GAS 
of A, W and their generalized inverses we have the following Jordan canonical 
forms, 


WA = Q 


VViiAil U 

0 W^ 22 A 22 


Q 


(WA)" = Q 


0 


0 


[(WA)“+']+=q[ "iQ 


-1 


^ J -6 Q 0 J ^ ■ 

Simple computing of W+ |(WA)^ [(WA)^^"*"^] ^ (WA)^| proves the theorem, 


W+ |(WA)'= [(WA)2'=+i]+ (WA)'=| = 

W-i 0] [(WnAn)'= O] [(WnAn)-^'^-! 0 

0 0 
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0 0 


0 


0 


■(WnAn)'= O' 

0 0 


■Wri'(WnAn)'=(WnAn)-2'=-i(WnAn)'= 0 

0 0 


Q ' = 


0 

■Wri^(WiiAii)-i 0 

0 0 

(WiiAiiWii)-i 0 

0 0 


= 


Q ' = 


Q ^ = Arf^w 


Using (1^ . an entry a^’^ of the W-weighted Drazin inverse A^^w can be ob¬ 
tained as follows 

m m m . 

-r "USE ».‘i’ (»;?"■’) (») 

s=l /=l 
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for all z = 1, m, j = 1, n. Denote by w^. the tth row of W* =: W = {wij) ^ 
]fjmxra for alH = 1, m. It follows from ^ . and (|T5|) that 


E rdetaWW*),.(wr.) 


E 


(k) + 

n;, w7 


'"ti 

1=1 1=1 


(fc) ae-fri,nb} 


E l(ww*) 

OL^Ir^ , n 


n{i} 


E rdetJ(WW*), (wt. 


E I(WW*) 

Ot.€lr^ , n 


(31) 


where ri = rankW. Similarly, denote by Vi. the tth row of =: 

V = {vij) S for all t = 1, m. It follows from E 

(fT5|) that 


771 _L 

S = 1 


(k) 


E rdett((v2'=+i(V2'=+i)*)^ 


j2 I (V2'=+i)*) 


OcGir 


E rdethfv2fc+i(V2fc+i)*) (v,.)) 

le/r, m{i} / t. / 

E |(V2fc+l (V2'=+l)*) “I 

aG/r, m 


(32) 


where r = rankW^+^ = rankW^. Using (IHTl) and (15^ in (1501) we obtain the 
following determinantal representation of Ad,w > 


d,W 

% = 


E E rdeu((v2'=+i(V2'=+i)*) (v.,))“ E rdeth(WW*)^.. (wO) 

t=laelr,m{t} ) t. I ae/,,^ „{j} ^ ^ 


E |(v2fc+i (V2fe+i)*) “I E l(ww*) 


aG/r, 


a€lr-i, 


(33) 


Thus we have proved the following theorem. 

Theorem 3.4 Let A e andW € withk = max{ Jnd( AW), Jnd(WA)} 

and r = rank(AW)^“*“^ = rank(AW)^. Then the W-weighted Drazin inverse of 
A with respeet to W possesses the determinantal representation i3S\} . where 
V = AW, V = v'=(V2'=+i)*, and W = V'^W*. 

Similarly we have the following theorem. 
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Theorem 3.5 Let A G andW G withk = max{Jn(i(AW),/n(i(WA)} 

and r = rank(WA)^“*“^ = rank(WA)^. Then the W-weighted Drazin inverse of 
A with respect to W possesses the following determinantal representation, 


d.W 

% = 


t E cdet, ((W*W) ^ (w.t)) ^ E cdett(((U2'=+i)*U2'=+i) 

*=l/3eJri.m{d /3GJr,„{t} 


E (w*w)^ E ((U 2 '=+i)*u 2 '=+i) ^ 

P&Jri.m ySeJr, n 


/3 

/3 


(34) 


w/iere U = WA, U = (U2'=+i)*U'=, and W = W*U'=. 

Proof. Using (l29)l . an entry of the W-weighted Drazin inverse A^^w can 
be obtained as follows 

= ttt (“!f+■') blf (35) 

s=l /=1 


for all i = l,m, j = l,n. Denote by w.t the tth column of W*U^ =: W = 
(wij) G for all t = l,n. It follows from E'''^*s'“i? = ^ t and (fTH) that 




E cdet, (W*W), (w*J ^ 

/3eJri,m{^} (k) 


S =1 


.=1 E (w*w) ^ 

^GJri, m 


uW = 


E cdeU((W*W).Ew.t))^ 

/3eJri, m{i} 


E 


(W*W) 


, (36) 


where ri = rankW. Similarly, denote by u ^ the jth column of (U^^+i)*!!^ =: 
U = {uij) G for all j = l,n. It follows from E = u.j and 
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(dH) that 


n , 

■>)»«’ = 


E 

1^1 


E cdett(((U2'=+i)*U2'=+i) 

PeJr,n{t} ^ ^ ^ 


Y, ((U 2 fe+i)*u 2 '=+i) ^ 

;SeJr.n 


(fc) 

•<■ = 


E cdet* (((U2'=+i)*U2'=+i) (u,)) 

3eJr,„{t} / .t / 


E ((U2fc+l)*U2fe+l) ^ 


(37) 


/3eJr, n 

where r = rank(AW)^+^ = rank(AW)^. Using the equations (l37ll and (l36ll in 
(|35|l . we obtain (l34|) . ■ 


3.3 Determinantal representations of the W-weighted 
Drazin inverse in some special case 

In this subsection we consider the determinantal representation of the W-weighted 
Drazin inverse of A € with respect to W € in a special case, when 

AW = V = {vij) G and WA = U = (uij) G are Hermitian. Then 

for the determinantal representations their Drazin inverse we can use m and 
m- For Hermitian matrix, we apply the method, which consists of theorem on 
the limit representation of the Drazin inverse, lemmas on rank of matrices and 
on characteristic polynomial. This method was used at first in |35] and then 
in [laiMlIMlISI]- By analogy to the complex case [35] we have the following 
limit representations of the W-weighted Drazin inverse, 

Ad,w = lim (AI„ + (AW)'=+2) (AW)'= A (38) 

and 

Adw= limA(WA)'= (AI„ -f (WA)'=+2)"^ (39) 

where A G IR.+ , and R+ is a set of the real positive numbers. 

Denote by and the jth column and the *th row of respectively. 
Denote by V'= := (AW)'=A e and W = WAW G 

Lemma 3.1 If AW = V = (vij) G with Ind'V = k, then 

rank (V'=+2) . < rank (V'=+2) . (40) 

Proof. We have V'=+2 = V'=W. Let V,s(-Wjs) G (s ^ i), be a 

matrix with —Wjs in the (i,s) entry, 1 in all diagonal entries, and 0 in oth¬ 
ers. The elementary matrix (—Wjs), {s ^ i), is a matrix of an elementary 
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transformation. It follows that 


E«!? 


{v*+2) _(v<«).np,.(-%.) 


s^i 


Wsl 




_(fe) 


-(fc) - -(fc) 

2 ^ VmsWsl ■ . ■ 

i—th 


-(fe) - 

Vis Wsm 

S#j 

—(fc) — 

2 ^ VmsWsm 
s#j 


We have the next factorization of the obtained matrix. 


/ ^ ii(fc) 


s/i 

-(fc) - 

2 ^ vWsWsl 


'C. 


If 

(fc) 


mj 

i—th 


E \ 

— (^) - 

2 ^ VmsWsm 


/ -(fc) -(fc) 

' v)-. v)c^ 


11 


12 


V. 


(fc) .-.(fc) 


21 


22 


V V. 


(fc) -(fc) 

ml Vm2 


Denote W := 


Wll 


V Wnl 


-(fc) 
Kn 
-(fc) 
V2n 

-(fc) 
Vmii ) 


0 

1 

0 

i—th 


( Wii 
0 

V Wnl 

Wlm 

0 

Wn.rn. 


0 

1 

0 

i—th 


Wir, 


j - th. 


j — th. The matrix W is obtained 


/ 


from W = WAW by replacing all entries of the jth row and the *th column 
with zeroes except for 1 in the {i,j) entry. Since elementary transformations of a 
matrix do not change a rank, then rank , rank W |. 

It is obvious that 

rankV^ = rank(AW)^A > rank(AW)^+^, 
rank W > rank WAW > rank (AW)'=+2. 


From this the inequality (HOI) follows immediately. ■ 

The next lemma is proved similarly. 

Lemma 3.2 //WA = U = (u^) G with IndXJ = k, then 

rank(U'=+^)^ < rank (U'=+2) , 

where U'= := A(WA)'= G 

Analogues of the characteristic polynomial are considered in the following two 
lemmas. 


14 























Lemma 3.3 //AW = V = {vij) G jj^x™ Hermitian with Ind'V = k and 
t G R, then 

cdet, (AI™ + V'=+2) ^ = cp^A"-! + + ... + (41) 

where Cn^'^ = cdeti andc^s^^ = ^ cdeti ^(V^+^) . ^ 

0^Js, n{*} 

for all s = l^n — i,j = 1, n. 

Proof. Consider the Hermitian matrix (tl + V^+^) ^ G Taking 

into account Theorem l2.8l we obtain 

det (AI + V'=+") ^ (v,^ = diA"-i + d2A"-2 + ■ ■ • + d„, (42) 

where ds = \ is the sum of all principal minors of order s that 

/3eJ..n{i} _ 

contain the i-th column for all s = 1, n — 1 and dn = det (V^+^). Consequently 

= where is the Zth column- 


vector of = (AW)*^A and WAW = W = (wu) for all I = l,n. Taking into 
account Theorem 12.11 we obtain on the one hand 

det (AI -b V'=+2) ^ = cdet, (AI -b V'=+2) , = 

= X] cdeti (AI -b V'=+2) j = X^cdeH (AI -b V'=+2) . ■ wu 

(43) 

On the other hand having changed the order of summation, we get for all s = 

1, n — 1 


we have = 


/ ^ ^11 ^Ll 

I 

V -{k) - 

Z2^21 


V -{k) - 

> v\,Wu 


ds= E det (V'=+2) I = E cdet. (V'=+2) ^ = 

/3eJs,n{d /3eA,n{i} 

E Ecdec((V'=+2) 7v«u;h))^ = 

E E cdeti ((V'=+2) . (v^f^)) ^ -Wii. 


By substituting (H51) and (H41) in (H^ . and equating factors at wu when I = j, 
we obtain the equality (HD) . ■ 

By analogy can be proved the following lemma. 

Lemma 3.4 //WA = U = (uy) G is Hermitian with IndU = k and 

t G M, then 

rdetj-(AI -b U'=+2)j , (uf ^) = -b ... -b 
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wherern^^ = rdetj(U^+^)j, andri^^^ = J 2 rdetj ^(U^+^)j. (u-^^)^ 


a€lB,n{j} 


for all s = l,n — 1 and i,j = 1 , n. 


Theorem 3.6 If A £ W S and AW = V = (w^-) € is 

Hermitian with k = max{Jn(i(AW), IndfWA)} and rank(AW)^+^ = rank(AW)^ 

r, then the W-weighted Drazin inverse Ad = (a’lf) e with respect to 

W possess the following determinantal representations: 


^E^^^cdet.((AW)^f (v'f)) ^ 
E |(aw)"+2 ^ 

h&Jr-, m ' 


(45) 


where is the jth column of'V^ = (AW)*A for all j = l,m. 

Proof. The matrix (AI^ + (AW)^+^) ^ £ jjmxm full-rank Hermitian 
matrix. Taking into account Theorem 12.41 it has an inverse, which we represent 
as a left inverse matrix 


(AI^ + (AW)'=+") ' = 


( ill L 21 . . . Lml \ 

Li 2 L22 ■ • • Lm 2 


det(AI™ + (AW)'=+2) 


\ I^lm 1^2771 • • ■ I^rnTTi j 

where is a left ij-th cofactor of a matrix Aim + (AW)^+^. Then we have 
(Aim + (AW)'=+2)"^ (AW)'=A = 

/ ^ . _(/c) ^ . _(fc) ^ . _(fc) \ 

/ . LslVg^ 2^ ^slVg2 • • ■ LslVsn 


det(AI„-|-(AW)''+2) 


s=l 


s=l 


s=l 


E l.2€;^ 




E l.2€^ 


(k) 


s=l 


(fc) 


E LsmVsl 

\ s=l 


E Ls 




S = 1 


... Ls2Vsn 

s=l 


^ r -(fe) 
• • • ^ J-^sm'^sn 

s=l 


By (1551) and using the definition of a left cofactor, we obtain 


Ad,w 


( 


lim 

CK^O 


V 


cdeti(AI™-|-(AW)''+2) 

det(AI,„-|-(AW)'=+2) 

cdet„(AI„,-|-(AW)''+^) 

det(AI™-|-(AW)'=+2) 


By Theorem 12.81 we have 


cdeti(AI„-|-(AW)''+=) ^ (v(^>) ^ 

det(Allir+(AWpT2) 

cdet„(AIm-|-(AW)*^+^) 

det(AI™-|-(AW)'=+'2) / 

(46) 


det (Aim + (AW)'=+2) = A’" + diA™-^ + dzA’"”" + .. • + dm, 
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where (is= (Aim + (AW)^+^) ^ is a sum of principal minors of (AW)^+^ 


3eJ. 


of order s for all s = l,77i — 1 and dm = det(AW)^+^. Since rank(AW)^+^ = 
rank(AW)^+^ = rank(AW)^ = r, then dm = dm-i = ... = dr+i = 0. It follows 
that det (AI™ + (AWf+2) = A"* + diA™-i + daA'"-^ + ... + Using 

m we have 


cdet^ (Aim + (AW) 


k-\-2 




(ij) yn-2 


niv) 


for i = l,m and j = l,n, where ^ cdeti ^ 


deJs,m{i} 

all s = 1, m — 1 and Cm'^ = cdeti(AW)^^^ 


m 


We shall prove that ^ = 0, 
when k > r+1 for i = l,m and j = 1, n. By Lemma lOl 

then the matrix more r right-linearly independent 

columns. Consider (^^0) ^ Js,m{i}- It is a principal 

submatrix of ^)) o^'^er s > r + 1. Deleting both its Uth row 

and column, we obtain a principal submatrix of order s — 1 of (AW)^+^. We 
denote it by M. The following cases are possible. 

• Let s = r + 1 and det M ^ 0. In this case all columns of M are right- 

linearly independent. The addition of all of them on one coordinate to 
columns of d ^®®P® their right-linear independence. 

Hence, they are basis in a matrix 

umn is the right linear combination of its basis columns. From this by 
Theorem 12.31 we get cdet^ ^(AW)^+^ (^3^0) /® ~ ^ <^s,n{*} 

and s = r -I- 1. 

• If s = r-|-1 and det M = 0, than p, {p < s), columns are basis in M and in 

(^d^)) d' Thenbv Theorem lOl cdet, ^(AW)^+^ (^A^)) d ~ 

0 as well. 


• If s > r -I- 1, then detM = 0 and p, {p < r), columns are basis in the 
both matrices M and A Theorem 12.31 we 

also have cdeti (^A^)) d 

Thus in all cases we have cdet^ ~ when P € Js,m{'i} 

and r -|- 1 < s < TO. From here if r + 1 < s < to, then 

cto) = ^ cdet. ((AW )'=+2 (v«)) ^ = 0 , 

/3eJ„,m{i} 
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and Cm ^ = cdeti ~ ^ ~ 

Hence, cdeti (AI + (AW)'=+2) ^ . .+ci"^'> 

for i = 1, TO and j = 1, n. By substituting these values in the matrix from (j46l) . 
we obtain 




Ad,rv = lim 


A-s-O 


"+...+c; 


( 11 )) 


(In) , m-1 


(ln),m-r 


A™+di A™-i + ...+dr'^"*“'' 


__+...+c;_ 

A’^+di A^“i + ...+dr 




+...+c: 


(ml)Am-r 


\ ii_ 

^ A^-j-fiiA^ ^-|-...-|-d7’A^ ^ 

/ 


(mn) , 




A^-j-d^A^ ^-j-...-j-fiT'A”^ ^ 


V 


dr 
,(n.l) 


dr 

.Cmn) 

_ 

dr 


Here ^ 

PeJr,m{i} 


cdetj ((A'=+i) ^ ^ and dr = J2 

^ ■ \ // BeJrr. 


(A^+i) 


Thus, we have obtained the determinantal representation of Ad,w by (1451) . 
By analogy can be proved the following theorem. 


Theorem 3.7 If A £ W e IHI"^^’”, and WA = U = (m„-) S is 

Hermitian with k = max{Ind{AW), IndfWA)} and rank(WA)^+^ = rank(WA)* 
r, then the W-weighted Drazin inverse Ad^w = ^ ijmxn respect to 

W possess the following determinantal representations: 


d,W 


E rdet, ((WA)^^+ 2 (uf))) “ 
eAnfH ^ ’ 

E |(WA)'=+"g 

ae.fr, n 


(47) 


where 


is the ith row o/U^ = A(WA)^ for all i= l,n. 


4 An example 


In this section, we give an example to illustrate our results. Let us consider the 
matrices 


A = 




/O i 
fell 
1 0 0 
Vl -k -jJ 


/ k 0 i 0 \ 
W = -j yfc 0 1 

\0 10 -kj 


Then 


V = AW = 


/ -k 

-1 -d 
k 

\-i + k 


-J 0 
i + k j 
0 i 
1-j i 


* \ 

1 + j 

0 

i — kj 


(* i 

U = WA = 0 k 

\0 0 
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and rankW = 3, rankV = 3, rankV^ = rankV^ = 2, rankU^ = rankU = 2. 
Therefore, IndV = 2, IndV = 1, and k = max{/nd(AW),/nd(WA)} = 2. 
It’s evident that obtaining the W-weighted Drazin inverse of A with respect to 
W by using the matrix U by p4l) is more convenient. We have 

/-I i + k 0\ fi 2+ 3j 0 \ 

= 0 -1 1 u® = 0 k 

\ 0 0 0 / \0 0 0 / 


/ -i 0 0\ 

(U5)*= 2-3j -k ,(U5)*U5 

\ 0 0 0 / 


/ 1 -2i- 3k O' 

2i + 3fc 14 

Vo 0 0 


u = (u^yv^ 


w*w = 


(^ 

—i 

j 

\-J 


( * 

1 + j 


f-k 

0 

j 

-k 

o\ 

1 

-2 + 3, 

—i + 6k 


—i 

0 

0 

V 0 

0 

n 

0 / 

u 

1 

k) 


i 

-j 

j ] 


f-k 

l-2j 

o\ 

2 

0 

-2k 

,w = W*u 2 = 

0 

i + k 

0 

0 

1 

0 

i 

1 +J 

0 

2k 

0 

2 ^ 



-1 

0 / 


Since by (l34l) 


d,W 


E E cdeti ((W*W) 1 (w.O) ^ E cdett(({V^)*vA (u, 

t=i/3e/3,4{i} /3eJ2,3{t} ^ 


E 

0^J3, 4 


(W*W) 


E 

/ 3 G J2, 3 


((u^rus)^ 


where 


cdeti {(W*W) , (w 1 )) ^ = 

/3e/3,4{i} 


k 

i 

-n ( 

cdeti 0 

2 

0 + cdeti 


0 

1/ \ 


^ cdeti ((W*W) , (w, 2 )) I = 
/3e/3,4{i} 


i 

d \ 

k 

-j 

j\ 

2 - 

-2k + cdeti 

r 

1 

0 = 0 , 

2k 

1 ) 

VO 

0 

2 / 

-2+ 

E cdeti 

((W*W) 

1 (w. 3 )) p 


/3e/3,4{i} 


= 0 , 


E |(w*w) 

4 


= 2 , 
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and 


E cdeti(((U5)*U5)^(u.i))^ 

/3e ^2,3{i} 


cdeti 


i —2i — 3fc 
-2 + 3j 14 


■ cdeti 


i 0 
0 0 


= 


cdet2(((U5)*U5) ^(u.i)) ^ = 0, 

/ 3 e ^2,3(2} 

^ cdet 3 (((U 5 )*U^) (u.i))^ = 0 , ^ |(m*U^) 

/3GJ2.3{3} /3GJ2,3 


= 1 , 


then 


d,w (0 • t) + (—2j • 0) + (0 • 0) 

«„ =--= 0 . 


Continuing in the same way, we finally get, 


^d,W = 


/ 0 -i 0\ 

0 0 0 

-1 5i-2k 0 

V 0 0 oy 


(48) 


By (HI we can obtain 


^-i -3 + 2j 0\ (-i -5 0> 

(U®)^=(0 -k , (WA)^ = = 0 -k 


0 0 Oy 

We verify (H5)) by (fTTl) . Indeed, 

k 0 i 0 

WAd,w = \ —j k Q 1 

0 1 0 -fc; 


0 0 0 ; 


/ 0 -i O' 
0 0 0 
-1 5z-2fc 0 


/-* 

-5 

o\ 

= 0 

-k 


Vo 

0 

oy 


= (WA) 


D 


We also obtain the W-weighted Drazin inverse of A with respect to W by (l24l) . 
then we have 


/O 


Ad,w = A ((WA)^)) = 


— % 




-k 6 + 5i 0 
-1 5* + 5fc 0 

\-i 5* + 6fc oy 


(49) 


The W-weighted Drazin inverse in (1491) different from (l48l) . It can be explained 
that the Jordan normal form of WA is unique only up to the order of the Jordan 
blocks. We get their complete equality, if Ad^w from (|4^ be left-multiply by 
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the nonsingular matrix P which is the product of multiplication of the following 
elementary matrices, 


P = P2,4{-k) ■ P4.3(-l) • P3.4(-6) • P4.r(-j) 


/1 

0 

0 


0 

1 

0 

-k 

0 

0 

7 

-6 

\-J 

0 

-1 

1 / 


Note that we used Maple with the package CLIFFORD in the calculations. 
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